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Abstract. We have studied the high-frequency properties of the non-equilibrium 
electron gas in GaN samples subjected to electric and magnetic fields. Spectra of 
the complex tensor of the dynamical mobility have been calculated for THz frequency 
range. For the compensated GaN and low temperatures, in the intervals of electric 
fields of the few kV/ cm and magnetic fields of the few T the existence of the cyclotron 
and optical phonon transit-time resonances has been identified. We have shown that 
interplay of two resonances gives rise to specific spectra of THz transmission and 
absorption (or gain). We suggest that experimental investigation of these effects will 
facilitate elaboration of field controlled devices for THz optoelectronics. 



PACS numbers: 72.20.Ht, 72.20.Dp, 85.30.-z 
accepted for publication in Semicond. Sci. Technol. 



High-frequency response of GaN 



2 



1. Introduction 

Progress in the technology of group-Ill nitrides, discovery of unique material properties 
and perspectives of high-power, and high-frequency applications [1] have inspired studies 
of high-field transport regimes in these materials [2j |3j H], including the high-frequency 
phenomena [5j El [7]. Most of these studies were focused to very high electric fields 
(10. ..200 kV/cm), where the intervalley electron transfer and the Gunn effects are 
expected. Also there is a necessity of understanding of hot electron kinetics in the 
range of moderate electric fields (1..10 kV/cm), because of fundamental and practical 
importance of such a case. Indeed, it is believed (see recent review [E]) that in polar 
nitride materials subjected to a moderately high electric field, strong inelastic optical 
phonon scattering can provide the so-called streaming transport regime. The latter 
is characterized by considerable anisotropy of electron distribution in the momentum 
space, quasi-saturation of the current-field dependence, negative dynamic conductivity 
and other useful effects. 

In general, the following conditions are favorable to realize the streaming regime: 
(i) strong electron-optical phonon coupling, (ii) low temperatures which provide optical 
phonon emission as dominant scattering mechanism of hot electrons and (iii) relatively 
low electron concentrations, when e-e scattering is negligible. According to the analysis 
given in a number of the works [E] , these conditions can be met in the nitride materials, 
as well as in the nitride heterostructures [9j [10]. At the streaming regime, the electron 
motion is nearly periodic: in an appropriate range of the dc electric field, an electron 
accelerates quasiballistically until it reaches the optical phonon energy, then the electron 
loses its energy by emitting an optical phonon and starts the next period of acceleratio.fi 
Such a mechanism of the nearly periodic electron motion leads to electric field induced 
resonances in high-frequency conductivity [the so-called optical phonon transient time 
resonance (OPTTR)] [13]. The OPTTR provides an operating principle for electrically 
pumped semiconductor sub-THz and THz sources [13], that has been practically realized 
in n-InP crystals [T4] . 

Experimental observation of the OPTTR in the nitrides and its further device 
implementation is impeded by the deficit of electron transport data for the range of 
moderate electric fields, including high-frequency properties. 

Additional useful information on properties of the nitrides can be obtained by 
investigation of magnetic field effects on charge carriers, for example, investigation of 
the electron cyclotron resonance (CR). It is well known, that exploitation of the CR is 
a powerful spectroscopic method that allows to determine important parameters of the 
band structure and carrier dynamics (for a detail review see ref. [IS]). Monitoring the 
cyclotron resonance in an electric field, one can study intraband dynamic processes 
of nonequilibrium carriers. For example, the cyclotron resonance technique was 
successfully used to study hot electrons and holes in Ge [16], Si [17], InSb [18], GaAs [15] 
and other important semiconductor materials. As the result, the detailed features of 

% Note, this simplified picture can be traced back to papers by Shockley pj] and Price [12] . 
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hot carrier kinetics have been understood. These include: relaxation mechanisms and 
their parameters, character of carrier distribution in the momentum space, the dynamic 
conductivity, etc. Recently, the cyclotron resonance was observed in the nitrides at 
moderate and high magnetic fields (above 5 Tesla). Here the resonance became apparent 
in THz-frequency range pJ2, EQj EU E2] , because of small sub-picosecond relaxation times 
characteristic for the nitrides. 

In this paper by the use of the Monte-Carlo method, we study GaN bulk-like 
material in moderate electric and magnetic fields. We focus on high-frequency response 
of GaN under conditions, when both resonances, the CR and the OPTTR, become 
apparent simultaneously. We show that overlapping of the two resonances produces 
specific spectra of THz transmission and reflection, as well as characteristic polarization 
effects. 

2. Model of electron transport 

We consider a bulk GaN sample of the cubic modification with a given concentration of 
ionized impurity, iVj. To exclude quenching effect on the OPTTR by electron-electron 
scattering we suppose that the sample is compensated, thus the electron concentration 
n e < iVj. At electric fields of a moderate strength, all electrons remain in the T valley 
and can be characterized by the parabolic dispersion law with the effective mass, m*. 
The dc electric field, F, and the magnetic field, B, are assumed to be parallel. Let the 
both fields be directed along OZ axis, F || B || OZ. The magnetic field is treated as 
a classical one. The ac electric field is supposed to be harmonic, f^exp^—iut) , with u 
being the frequency. To find the small-signal response, the ac field is considered as a 
small perturbation, \f u \ <C |F|. 

To calculate electron transport characteristics and, particularly, frequency 
dependent dynamic mobility, we use the single-particle algorithm of the Monte-Carlo 
procedure [23], [25] . For the simulation of the electron transport we take into account the 
main three scattering mechanism: scattering by ionized impurities, acoustic phonons 
and polar-optical phonons. For the ionized impurity scattering, we exploit the "mixed" 
scattering model that unifies both Brooks-Herring and Conwell-Weisskopf approaches. 
This approach is more appropriate for the analysis of compensated materials. The 
formulae of electron scattering probabilities for these mechanisms can be found in 



In the limit of a classical magnetic field the scattering probabilities do not depend 
on this field. However, the magnetic field determines the Lorentz force, so that electron 
dynamics between sequential collisions is described by the equation : 



where v is the electron velocity, e is the elementary charge, c is the light velocity. To 
avoid cumbersome expressions for v(t) during free flight between collisions we present 



ref. [23J. 



m*v = — e F + f w exp(— iut) H — [v x B] 



(1) 
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here only velocity modulation by the ac field, f w : 




sm(u c At) + 
cos(o; c At) — 




(2) 



where to is an initial moment, At = t — to and u c = eB /cm* is the cyclotron frequency. 
We used Eqs. §B) and ([2]) to apply the Monte-Carlo procedure to the steady-state kinetics 
and high-frequency response of the electron gas subjected to both electric and magnetic 



First, we studied the steady-state distribution function of hot-electrons, the drift 
velocity, Vd(F), the average energy, E(F), and other stationary electron characteristics 
in wide range of applied dc electric fields, when the magnetic field, B is present. In FigJT] 
we show an example of calculations of Vd(F), E(F) and distribution function. These 
results are discussed below in more details. 

For the magnetic fields of a moderate strength (B < 10 T), we found that 
these characteristics practically are independent on the magnetic field. This is a 
consequence of two main factors: (i) for the parabolic electron dispersion and the parallel 
configuration of electric and magnetic fields, electron motion along fields direction and 
that in perpendicular plane are uncoupled; (ii) the magnetic field does not affect the 
energy balance of the electron system. The detail analysis of the steady-state transport 
characteristics in compensated GaN can be found in ref. [24] . 

In the presence of steady-state electric and magnetic fields, the high-frequency 
response can be characterized by the dynamical mobility, which is the tensor, (i^. 
The dynamic mobility tensor defines the ac electric current induced by the ac field, 
J w = euefi^. For the considered configuration of the fields, this tensor has five non- 
zero components: 



where only three components, ^^ x , ^ w ,x y and fi UjZZ , are linearly independent. Other 
components are \i^ m = fJ-ui,xx, Hu,xy — —^ui,yx- Note, at B = nondiagonal components 
are equal to zero. In next Section we present the analysis of spectra of all components 
of fj, u and their dependencies on the electric and magnetic fields. 

3. Calculation of the dynamic mobility 

Prior to discussion of results on the hot-electron dynamical mobility, £i w , we shall 
estimate ranges of the parameters favorable for observation of the OPTTR and the 
CR effects in GaN. 

For concentrations of the electrons and ionized impurities we set n e = 10 15 cm -3 
and Ni = 10 16 cm -3 , respectively. The lattice temperature is assumed to be 30 K. For 



fields [231. 




(3) 
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the GaN material parameters, such as effective mass, dielectric constants, deformation 
potential, etc, we used the database from ref. [26] The low field mobility was calculated to 
be ~ 5000 cm 2 /V s, which is well correlated with the measured mobility in GaN samples 
of low dislocation density [27] . For these parameters, the streaming transport regime and 
the OPTTR are the most pronounced in the range of the dc electric fields 1...10 kV/cm. 
As seen from Fig. [TJ (a), in this range of applied dc electric fields the drift velocity 
and the averaged energy show the typical quasisaturation behavior. As the results of 
streaming electron motion, the strongly asymmetric electron distribution function is 
formed in momentum space. In Figs. [I] (b) and (c), we present both longitudinal and 
transversal profiles of the distribution function for the case of developed streaming, 
F = 3 hVj cm. As seen, the transversal profile of the distribution function remains 
symmetric. While for the longitudinal one shows strong asymmetry: the electrons have 
momenta mainly directed along the applied force. The streaming distribution leads to 
the effect of OPTTR in the frequency range 0.2 — 2 THz [8]. 

Classical treatment of electron motion in the magnetic field requires the condition 
hu c <C E. For the above range of the electric field, the average hot electron energy 
in GaN is estimated to be E < 30meV. Then, the latter inequality gives limitations: 
u c <C 45 THz and B <C 50 T. In what follows, we restricted ourself to the analysis 
of moderate magnetic fields, 1...5T, which correspond to the cyclotron frequency 
w c /2tt « 0.1...0.7THz. 

3.1. General properties of fi^ tensor 

A representative example of the dynamic mobility of hot electrons under the developed 
streaming regime (F = 3 kV/ cm) is given in Fig. [2j The magnetic field is set to be 
B = 4.5 T. In this figure, three independent tensor components, \i w ^ xx ^ ^u),zz, Hw,xy, are 
shown as functions of the frequency, u. 

The 'longitudinal' (with respect to the applied fields) component, fiu, lXX , is 
practically independent on the magnetic field and demonstrates all signs of the OPTTR 
effect: an oscillation behavior of both Re\jj, UtZZ ] and lTn\pL UjZZ ], and a frequency 
"window", where the real part of the dynamic mobility is negative, Re\jJL UlZZ ] < 0. 
Minimum of Re\p UjZZ ] reaches the value ~ — 200 cm 2 /Vs at u/2tt = O.Q3THz. Notice, 
absolute values of -Ref/i^J including the case u — Y are much smaller than the low 
field mobility. This is a consequence of above mentioned quasi-saturation of the velocity 
- field characteristics. 

The 'transverse' component, fM^^x, and off-diagonal component, fJL u , X yi are 
dependent on both F, B fields and demonstrate features proper to the CR. The 
maximum of Re\jJL U)XX ] and zero of Im\p u>xx ] occur at u — > u c , u c /2tt w 0.63 THz for 
the used parameters. In the vicinity of u c , we find that Re\pL W)XX ] is considerably larger 
than |i?e[// w>22 ]|. This is because all electrons contribute to the CR, while exclusively 
streaming electrons contribute to the OPTTR effect. 

We found that spectral dependencies of two tensor components, n w ,xxi AV^y can 
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be approximately described by the simple Drude-Lorentz model that takes into account 
electron motion in a magnetic field. For this model, the components are 

er* 1 — iur* 



IK 



(DL) 



JPL) = _tl ^ (A) 



m* (1 - iur*) 2 + (u c r*) 2 
er* uj c t* 
m* (1 — iur*) 2 + (uj c t* 

where r* is an effective relaxation time. As seen from figures |2] (b) and (c), Eqs. (Tj0) 
give a good fitting to the Monte Carlo result for r* = 0.6 ps. Herewith, the CR is well 
pronounced, because u c r* m 2.4. 



3.2. Electric- field dependence of the cyclotron resonance 

For the analyzed configuration of F and B fields, the electric field does not affect the 
cyclotron motion of the electrons in the {x, y}-plane. Instead, this field modifies the 
electron distribution and alters relaxation processes, and, consequently, the effective 
relaxation time, r*. 

The results presented in Fig. [3] illustrate the electric field effect on the cyclotron 
resonance. One can see that, when the F field varies, the frequency of the resonance 
remains equal to u c , while the sharpness of the resonance degrades considerably with 
increase of the F field magnitude. For example, at B = 2.2 T (u c /2tt = 0.32THz) the 
cyclotron peak in Re\j/, u>xx ] practically disappears when F > 5kV/cm (see Fig. E^a)). 
This is, obviously, due to an intensification of scattering processes. For the larger 
magnetic field, B = 4.5 T, the CR remains observable for the range F = 1...5 kV/cm, 
as shown in Fig. [3](b). Therein we present also the F field dependence of the effective 
relaxation time, r*, found by fitting Eqs. H]to the Monte-Carlo results. It is interesting, 
that for the F field range, 1..5 kV/cm, when the streaming regime is formed, t*(F) 
dependence qualitatively correlates with the transit-time, r tr , during which a ballistically 
accelerated electron reaches the optical phonon energy, fujo, i.e., r tr = ^/2m*hu)o/ eF . 
Thus under very fast process of optical phonon emission, the width of the CR is 
determined rather by the dynamic characteristic, r tr . 



4. Coexistence and interplay of the CR and the OPTTR 

To measure microwave/optical effects in a magnetic field one uses, typically, the Faraday 
and/or Voigt configurations of an incident wave and a magnetic field. The CR can be 
observed for both configurations, while the OPTTR can be measured only for the Voigt 
configuration (if F and B fields are parallel). To study coexistence and interplay of both 
resonant effects we will focus on the Voigt configuration. 

For certainty, we set that the incident electromagnetic wave propagates along the 
OY axis and the high-frequency electric field, f w , lies in the {x, z}-plane. The incident 
wave is supposed to be linearly polarized. Then, for the wave polarization along the 
OX axis (f^HOX) the CR is observable, but the OPTTR is absent. To the contrary, for 



High-frequency response of GaN 



7 



the wave polarization along the OZ axis (f w ||OZ) the CR is missed, but the OPTTR 
appears. For other polarizations, both resonances coexist and affect each other. 

Interplay of the two resonances becomes apparent, for example, in transmission 
spectra. These spectra comprise two contributions from electron and lattice subsystems. 
In additions, the spectra vary with the sample thickness because of reflection and 
interference effects. In Fig. H^a), we show the transmission spectra of a GaN sample 
in the THz frequency range, where the CR and the OPTTR appear. The electron 
parameters and magnitudes of the F, B fields are the same as in Fig. [2j The sample 
thickness is set 10 fim, the dielectric permittivity is 8.9. For comparison, the contribution 
of the lattice subsystem is shown separately. In accordance with the above said, 
the transmission spectra depend on the polarization of the incident wave. For the 
polarization parallel to the F, B fields, the CR is absent, however due to the OPTTR 
the electron contribution leads to an increase in the transmission coefficient (curve 1). 
For the perpendicular polarization, the CR decreases the transmission (curve 3), while 
there is no the OPTTR contribution. The transmission spectra for an intermediate 
polarization are illustrated by curve 3. 

In Fig. H(b), we present frequency dependencies of the loss/gain coefficients which 
append the data in Fig. @|a). The loss/gain coefficient defines decrease/increase of the 
electromagnetic wave energy per unit time due to interaction with the electrons in the 
GaN sample. One can see that characteristic for the CR decrease in the transmission 
coefficient is due to cyclotron absorption, while sample 'bleaching' under OPTTR is 
relevant to a gain of the electromagnetic wave. 

The loss/gain effects and their dependencies on the wave polarization can be 
analyzed as follows. Let 9 be an angle between the electric vector f w and the OZ 
axis. Then, the time-averaged ac power density, that is received by electrons from the 
electromagnetic wave, reads: 

f.1 2 - (5) 

A negative value of P^ indicates that an energy is transferred from the electrons to the 
electromagnetic field and vice- versa, a positive P^ means that the energy of the wave is 
absorbed by the electron subsystem. 

In the frequency window, where i2e[/i W)Z J < 0, the function P w {6) can be negative 
in some interval of the polarization angles 9. The condition P w (9) = restricts the 
parameters 9 ,u, where the gain due to the OPTTR occurs. The actual space of these 
parameters can be illustrated in the {9, o;}-plane, as presented in Fig. [5] for different 
values of the B field. From this figure it is seen that for B = the gain occurs in the 
sufficiently wide angle sector, \9\ < 35°. With increasing B this angle sector becomes 
narrow. For example, at H = 4.5 T the gain is possible only at \9\ < 15°. Thus, 
polarization dependencies of the gain effect are strongly affected by the CR. 



-en P 



^,xxsin 2 (9) + ^,zzCos 2 (9) 
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5. summary 

In this paper, we have investigated bulk-like GaN subjected to electric and magnetic 
fields of moderate magnitudes, kV/cm and 1..5T, respectively. The parallel 
configuration of the electric and magnetic fields has been assumed. With the use of the 
Monte Carlo method, we have calculated steady state and high-frequency characteristics 
of the compensated GaN at low temperatures. Particularly, we have focused on 
the situation, when two resonances, the cyclotron resonance and the optical phonon 
transient-time resonance, become apparent simultaneously. We have calculated spectral 
dependencies of all components of the tensor of dynamical mobility in THz frequency 
range for different magnitudes of the electric and magnetic fields. 

We have studied the CR for hot electrons, particularly, the electric field dependence 
of the CR broadening. For the analyzed range of the electric fields, where the streaming 
transport regime and the OPTTR effect are realized, we have found that the effective 
relaxation time, t*(F), determinative the broadening, correlates with the transit-time, 
T tr , during which a ballistically accelerated electron reaches the energy of the optical 
phonon. That is, at very fast process of optical phonon emission, the CR broadening is 
determined rather by the dynamic characteristic, r tr . 

We have determined spectral and field dependencies of the OPTTR. Then, we have 
found that at electric fields of a few kV '/ cm and magnetic fields of a few T, well-developed 
cyclotron effect and optical phonon transient-time effects may coexist in the frequency 
range of 0.5. .. 1 THz: they can be observed simultaneously for the Voigt configuration of 
the incident THz wave and the fields. We have shown that interplay of two resonances 
gives rise to specific spectra of THz transmission and absorption (or gain), as well as 
characteristic polarization effects. 

It is important to make a few remarks on experimental observation of the discussed 
effects in the THz frequency range. The spectral measurements in the THz range 
can be performed by using recently developed detection methods for ultrashort electric 
signals and generation techniques for sub-picosecond transients of electric fields. The 
advanced ultrafast detection methods mainly relay on two techniques: photoconductive 
sampling [28] and electro-optic sampling [291 EQ] • By the use of femtosecond laser pulses 
with high repetition rates, these techniques allow for the direct measurement of electric 
transients with sub-picosecond time resolution. As for the generation techniques of 
sub-picosecond electric field transients, the most studied are THz pulses emitted by 
transient photocurrents in a photoconductive switch and THz emission generated by 
nonlinear optical processes (optical rectification, four- wave mixing, etc). Additionally, 
time-resolved THz polarization and ellipsometry methods are available. The techniques 
of generation and detection of sub-picosecond electric field transients were applied 
to GaN films to measure real and imaginary parts of the small signal conductivity 
under equilibrium conditions [SI El [7]. Apparently, the discussed THz techniques are 
generally compatible with high-field measurement set-ups for magnetotransport in GaN 
films/plates. Particularly, these techniques can be applied to study the CR and OPTTR 
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for the most interesting Voigt configuration of the fields. 

Concluding, the presented results show that low temperature investigation of GaN 
samples in moderate electric and magnetic fields can provide a useful information on 
hot electron dynamic and relaxation, as well as reveal novel high frequency effects. 
We suggest that experimental investigation of these effects would facilitate elaboration 
of field controlled devices for THz optoelectronics, including THz emitters, amplifiers, 
electro-optical modulators [31], etc. 
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Figure 1. Panel (a) shows the field dependencies of the drift velocity, Vd(F), and 
averaged energy, E(F). Panels (b) and (c) show the longitudinal and transversal 
profiles of electron distribution function at p± = and at p z = 0, respectively, at 
F = 3kV/cm. Characteristic velocity, Vq = 4 x 10 7 cm/s, optical phonon energy 
Huj = 92 meV, po — y/2m*hujo is the characteristic momentum. 
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<n/2it,THz 

Figure 2. (a), (b) and (c): Real (solid lines) and imaginary (dashed lines) parts of 
the components /x WjZ ^, /i w , x;c and (i u ,y X as functions of the frequency for F — 3 kV/cm 
and B = 4.5 T. Thin lines are fitting curves defined by Eqs.|U 
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ml2n, THz 

Figure 3. Real part of n u ,xx as function of the frequency for different F and B fields, 
(a): B = 2.2 T; (b): B = 4.5 T. Curves 1, 2, 3 are calculated at F = 1,3,5 fcF/cm, 
respectively. Vertical dashed lines depict positions of cyclotron frequencies. In the 
insert we compare t*(F) and r tr (F) dependencies (squires and solid line, respectively). 
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Figure 4. Transmission coefficient ( a) and ioss coefficient (b) of a GaN piate. 
Electron parameters and magnitudes of the F and B fields are the same as in figure [5] 
The plate thickness is lOfim. Curve 1: f\\OZ, curve 2: f||OX, curve 3: the angle 
between f and the OZ axis is 45°. Dotted line shows the lattice subsystem contribution. 
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Figure 5. Isolincs P u (0) = for different magnetic fields and F = 3kV/cm. Curves 
1, 2, 3 correspond to H — 0, 2.2, 4.5 T, respectively. 



